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AN EXACT TRANSFORMATION FROM GEOCENTRIC TO GEODETIC 


COORDINATES FOR NONZERO ALTITUDES 


David R . Hedgley , Jr . 
Flight Research Center 


INTRODUCTION 


In space position measurement systems , it is often necessary to transform geo- 
centric to geodetic coordinates . Consequently , it is important to have an exact 
closed-form solution for the transformation that is free of embedded singularities . 

An exact solution is particularly important for determining missile positions . All 
the previous work in this area falls into two categories: formulas that are mathemat- 
ically exact but involve equations that become unstable in the neighborhood of their 
singularities, and formulas that are approximations that may or may not have inher- 
ent singularities. References 1 and 2 are of the former type, and references 3 to 5 
are of the latter . 

Although some of these methods are adequate under some conditions, they have 
limitations under certain circumstances. This paper presents an exact closed-form 
solution for the transformation that is free of singularities. The solution permits 
the evaluation of any of the existing methods under any conditions and should serve 
as a primary standard . 

The appendix contains computer subroutines that implement the ideas that 
follow and examples of input and output data . 


SYMBOLS 


Parenthetical 
a (A) 
b (B) 
d 


symbols are computer identifiers for variables . 
major axis of ellipsoid 
minor axis of ellipsoid 

the distance from any point in space to the ellipsoid 



e (E) 
f (F) 
f(x,y,z) 

g 

H 

h (ALT) 
(ALTS) 

Po 

P2 


X 

x.y.z 

(X0,Y0,Z0) 

^ 2 ’^ 2 ’^2 
(X2,Y2,Z2) 

Y 

Z 



eccentricity 
flattening factor 

general function to be minimized 
restraint imposed on f 

function which is composite of functions f and g 
altitude normal to the ellipsoid 
altitude computed from distance formula 

any point in space 


the point on the ellipsoid that is the minimum distance from the 
point Pq 

percentage of error in d 

X 


percentage of error in d 


y 


percentage of error in d^ 

coordinate axis that intercepts the Greenwich meridian 
coordinates of any point on the ellipsoid 
coordinates of 


coordinates of P 2 


coordinate axis that is in a direction normal to the plane deter- 
mined by the intersection of X and Z 

coordinate axis that intercepts the poles 


2 


a (V) 


arbitrary real constant to be determined 


e 

(ALONG) 

X (ALONG) 
w 

cp (ALAT) 

CO 


error in a 
longitude east 

longitude west 

geodetic latitude 
geocentric latitude 


ANALYTICAL FORMULATION 


Let a model of the earth be an ellipsoid given by the following equation (ref. 6): 

2 2 ,2 

£_ + ^ + = 1 ( 1 ) 
a a b 

where a> b. In addition, let (Xp.yQ.z^) be any point above the ellipsoid, as 
shown in the sketch below . 
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The line perpendicular to the surface of the ellipsoid from Pq is the altitude of 
Pq above the ellipsoid and hence is the shortest distance from to that surface 
(eq. (D). 

Therefore, by minimizing the distance from p^ to the surface of the ellipsoid, it 

is possible to acquire the coordinates (^ 2 , 3 ^ 2 ’ point on the surface 

that cause the distance defined by 

= ^ < 2 ) 

to be minimum . 

The Lagrange multiplier method (ref. 7) is used to implement this minimization 
where 

2 2 2 
f(x,y,z) = (x - + (y - y^) + “ ^o) 


and 


2 2 2 
/ , X ^ y ^ z 

gix,y,z) = — 9 + "V + 


2 2 ,2 
a a b 


Then 


( 2 


, - 2 2 \ 


It partial derivatives are taken in turn and each is equated to zero, the following 
equations result: 


9H „/ \ «2jc „ ,(^ / 2'\ 

^ = 2 (o: - Xq) - = 0 or X = Xq/{1 - a/a ) 

a 


If = 2 (y - y„) - ^ = 0 or y = y„/(l - a/a^) 


(3) 


(4) 


(5) 


2 2 2 
dH _ _ X y z 

da ~ 2 2.2 

a a b 


+ 1 = 0 


( 6 ) 


4 


Equations (3) to (5) can be substituted into equation (6) as follows: 


■|jCq/(i - a/a^)j /a^ - (l - a/a^)J /a^ - [zq/(i “ + 1 = 0 


or 


x.^(b^ - 2a + a^/b^) - - 2a + a^/b^) - z - 2a +a^/a^) + (a^ - 2a - a^/a^)(b^ - 2a + a^/b^) = 0 


"O'" J'Q 

Simplifying and collecting , 


{l,aW - 2(i/»2 * wy * (4 * a^/o' » 

+ + a^b^ - 


2 2 2 „ 
Z0« = 0 


(7) 


This is a quartic equation in a and can be solved in closed form (ref. 8) . The 
proper solution can be determined if the solution is restricted to the real zeros of 
equation (7) and that zero of equation (7) is chosen that causes d to be smallest . 

The appropriate solution of this quartic has no singularities . That is , a does not 
approach plus or minus infinity , because that would imply from equations (3) to (5) 
that x^, y„, and equal zero. Equations (3) to (5) do not become unstable, 

Z Z Z 2 2 

because the statement a = a or a = b implies that x^, y^t or approaches 
infinity. Since altitude is nonzero, a does not equal zero. 

Let Oq be that zero that satisfies the quartic equation. Then, from equations (3) 
to (5) , 


^2 = ■ "o) 

^2 = “ "o) 

1/2 

'■s = [(^2 - ^ (>’2 ■ '’o)^ - ^o)d 

cp = sin''(Zo-22y'>s 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 
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Although the foregoing formulas are mathematically exact, for computational reasons 

it is usually advantageous to normalize all distances , because this minimizes errors 

due to rounding off. However , errors due to rounding off are still possible. To 

circumvent this problem for altitudes not near zero, consider the following error 

analysis with respect to the computation of h , X , and (p . 

s 


Let e be a small error introduced in the computation of the zero of the quartic . 
In addition, let ~ ^0 ~ ^2 ~ 
by Op only, because a, b, oCq, y^, and are all exact. It is easy to show that the 
percentage of error, t^, of is as follows: 


o“o 




Then the error in 


d is introduced 
z 


= -1 + (1 + e/a^)[(b2 - - e)_ 


From equation (10) , h‘ 
is small, (b^ - 


t ,t 

X’ y 


ela. 


2 2 I I 

- Oq > b . Therefore b - > a 

- Oq - ^ 1. Therefore ^ 


0 


> U 


, and since s 


Similarly , 


Since equation (12) involves the ratio of two numbers with virtually equal 
percentages of error , latitude is undisturbed by a small error in the computation of 
Up due to rounding off. 

As noted above, this analysis is valid only if the altitude does not approach zero; 
otherwise approaches zero and the value |g/“q| niay approach infinity. The 

instability of |e/aQ| suggests that equation (12) may become unstable under these 
conditions . 


This final problem can be resolved. If P 2 (:xt 2 ,y2 > •^2^ point on the 

ellipsoid that intersects the normal drawn from (Xg,yQ,z^) , the following equa- 
tion (ref. 6) may be used: 


cp = tan 


tan (o/(l - e^) 


(14) 


where tan 00 = ^^^{^2 ^ ^2^) 


1/2 


From these relationships , <p may be determined 


2 ■ ^2 

immediately . The error , if any , in tan co introduced by the computational round-off 
error in is inconsequential, because it can easily be shown that the percentage of 

error in ^ 2 , y^, and Z 2 is virtually zero. In fact, the percentage of error in z^ and 

x^ or y 2 is given by ^b^ - - e) - 1 and - 1, 

respectively. Since e is small for any value of a^, the quantities 
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(b^ - a^/(^b^ - Oq - and (a^ " ~ “o ~ ®) 1 under any condi- 

tions. Thus, equation (14) is preferable to equation (12) for determining geodetic 
latitude for any altitude . 


However, even the smallest error in causes more than negligible error in 

the computation of h (eq. (11)) , because the error is amplified by the scaling 

factor used to restore the true value of h after computations are completed . This 
analysis suggests that it is preferable to use the following exact formulas from 
reference 6 to compute h instead of equation (11): 


Xq = l' 


a/(l - 


e sin cp ^ 


cos A, 


a /2 1 

r / 22 1 

Q = |a/a - e sin cp / + hjcos cp sin X 

= [a(l - e'^)/(l - e' sin^ , h] si 


Sin cp 


2 2 

where f- {a - b)/ a and e = 2f - f . Since latitude and longitude are known, h 
can be computed directly and the proper equation can be chosen to compute h to 
avoid division by values close to zero . 

In this way the geodetic coordinates cp , A. , and h are determined exactly both 

mathematically and computationally , Therefore , the computational accuracy of the 
method is the accuracy of the computer used . 


CONCLUDING REMARKS 


A mathematically exact method for computing geodetic coordinates from geocen- 
tric coordinates is derived . The computational accuracy achieved by using the 
method is as accurate as the computer used , The transformation provides a primary 
standard and makes possible the evaluation of any of the existing methods . 


Flight Research Center 

National Aeronautics and Space Administration 
Edwards, Calif., November 4, 1975 
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APPENDIX 


COMPUTER SUBROUTINES AND DATA SAMPLES 


The following computer subroutines implement the theory presented in the text . 
The subroutines were written in FORTRAN IV . 

The comments in the listings should be sufficient for their comprehension and 
modification . A sample input and output listing is provided to facilitate the verifi- 
cation of the correct FORTRAN code . 
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APPENDIX - Continued 


Subroutine GEOD 


SliBROUTINE GEOO (ELEf ELEH , ELES ♦ ^ ZI fAZI , A ZIS » PH I , FHIM , P HI Si 1 

»ALH,ALHM»ALMS,RAN,H*AUAT tALONG t AL T) 2 

C ^ 

C ^ 

5 C THIS SUBROUTINE COMPUTES THE ALTITUCE, GEODETIC LATITUDE AND 5 

C LONGITUDE GIVEN THE RANGE, A7IHUTH AND ELEVATION OF A FCINT 6 

C WITH RESPECT TO A RADAR SITE. Z 

C ® 

C ° 

10 C ELE,ELEM,ELES ARE THE CEGREES, MINUTES AND SECONDS RESPECTIVELY OF 10 

C ELEVATION OF THE TARGET WITH RESPECT TO A RADAR SITE. 11 

C 12 

C AZI,AZIM,AZIS ARE THE DEGREES, MINUTES AND SECONDS RESPECTIVELY OF IZ 

C AZIMUTH OF THE TARGET WITH RESPECT TO A RADAR SITE. 14 

15 C 15 

C PHI, PHIM, PHIS ARE THE CEGREES, MINUTES AND SECONDS RESPECTIVELY 16 

C OF THE GEODETIC LATITUDE OF THE RADAR SITE. 17 

C 1^ 

C ALM, ALHM, ALMS ARE THE CEGREES, MINUTES AND SECONDS RESPECTIVELY IE 

20 C OF THE LONGITUDE OF THE RADAR SITE. 20 

C 21 

C RAN IS THE RANGE OF THE TARGET WITH RESPECT TO THE RACAR SITE. ??. 

C H IS THE ALTITUDE OF THE RADAR SITE ABOVE THE REFERE^CE ELLIPSOID. Z3 

C 

25 C ALAT IS THE COMPUTED GEODETIC LATITUDE. 25 

C ALONG IS THE COMPUTED LONGITUCE, ^6 

C ALT IS THE COMPUTED ALTITUDE. 27 

C 28 

DIMENSION R (12 ) .ROOT (1 ?> 29 

30 C -0 

C DEFINE MAJOR AND MINCR AVIS OF ELtlPSCIO IN FEET, 31 

C ^2 

A=209Z583?. ^3 

R=20654892. 34 

35 F=fA-R)/A 35 

CE=^^F-F»*2 36 

C 37 

C NORMALIZE IN UMTS OF A. 38 

C 39 

40 7P=A 40 

A1=A 41 

H=H/7P 4? 

RAN=RAN/ZP 43 

A=A/ZP 44 

t.5 R = B/ZF 45 

GK=3.141S92653€ 46 

CCN=GK/180. 47 

C 48 

C convert angles FROM DEGREES TO RADIANS. 49 

50 C 50 

ELE=ELE* (ELEM/60M-(ELES/3600) 51 

PHI=PHI+(PHIM/60) PHIS/ 3600) 52 

AZr=AZH* (AZIM/60) + (AZI S/360 0) 53 

ALM=ALM+ {ALMM/60)f( ALMS/3600) 54 

55 ALM=360-ALM 55 

PHI=PHI*C0N 56 

ALM=ALM*CON 57 


9 



APPENDIX - Continued 


AZI=A7I»C0N 58 

ELE=ELE*CON 59 

6D C 60 

c ei 

C CCMFUTE THE TRANSFOF^'ATI CN CF TARGET COORHIKATES, E? 

C FROP FCLAR CCCROINATES TC CARTESIAN COORCINATES, 63 

C CENTERED AT TRACKING RADAR. 64 

65 c es 

XO=RAN»SIN(AZI)*COS ( ELE) 66 

70=RAN*SIN(ELE> 67 

Y0=RAN»C0S(A7I) *C CSC ELE) 68 

C o9 

70 C COMPUTE THE TRANSFORMATION OF THE RADAR SITE GIVEN 70 

C IN GEODETIC POLAR CCCROINATES TO CARTESIAN CCORDINATES 71 

C WITH RESPECT TO THE GEOCENTRIC CENTER OF THE EARTH. 72 

C 73 

EX = l/( i-CE*( SIN (PHI) ♦♦2)> 

75 EX=ABS(EX) 75 

EX=EX**. 5 76 

EX=EX*A 77 

EO= (EX + H)*COS (FHI)»CCS (ALM) 78 

FO=(EX + H)»COS (PHI )^SIN( ALM) 7 9 

80 GO= (EX»(1-CE)+H> ^SIN(PHT) 80 

C 8 1 

C CCMPUTES THE DIRECTION COSINES OF THE AXIS OF THE RADAR SITE. 82 

C "3 

U1 = -SIN (ALM) 84 

85 U2=C0S(ALM) 85 

U3=0 86 

Vi=-CCS(ALM)*SIN(PNI ) 87 

V?=-SIN( ALM)*SIN(PMn 88 

V3=COS(PHI) 89 

90 W1 = C0S(ALM)*CCS(FHI) 90 

W2=SIN(ALM) ♦COSCFHT) 91 

W3=SIN (PHI) 92 

C 93 

C CCMFUTE the TPANSFORMATICK CF 7 FE FCINT CR TARGET 94 

95 C FROM CARTESIAN COCRCINATES CENTERED AT TRACKING 95 

C RADAR TO GEOCENTRIC CARTESIAN COCPDINATES. 96 

C 97 

X5=XO^UlfYO^V1^70*Wl+EO 98 

Y5=X0^f2+Y0* V2470*M2+F0 99 

100 75=XO»U3+YO»V3*7C*W3+GO 100 

XO=X5 101 

YC=Y5 102 

7C=75 103 

C 104 

1G5 C CONSTRUCT THE COEFFICIENTS CF THE CUARTIC ECUATIOM 105 

C R(5) * X^»4 + R(4> ♦ X»*3 ^ R(3) » y**Z ♦ 106 

C R(2)»XfR(l)=0 107 

C 108 

R(5)=l/( (A»R) »»2) 109 

110 R (4 ) = -2» ( (1/ (A»*2)M- (1/ (9»»2) ) ) 110 

R (3) =4+( (R/ A) ^»2) ♦(( A/D) ♦♦?)-( (XO/B) *♦21- ( (YO/e )**?)- C (70/AI^^»2) 111 

R (2)=2^(X0^^2 4Y0»»2 + 70*»2)-2» (A»*24D»^2) 112 

R (1 )= (A»D)**2-(X0»D) »^2-( Y0^P)^*2*(70»A)»*2 113 

C 114 
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APPENDIX - Continued 


115 


120 


125 


130 


135 


It+O 


145 


150 


1 55 


leo 


165 


170 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


CALL QUART TC SOLVE THE QUARTIC EQUATION FOR ALL 115 

REAL ZEROS WHEREt 116 

R = ARRAY OF COEFFICIENTS 117 

ROOT = ARRAY OF REAL 2ERCE 118 

NI = NLMBER OF REAL ZERCS 119 

120 

CALL QUART(R,ROOT,NI ) 121 

1^2 

OETERHINE THE COORDINATES ON THE ELLIPSOID THAT 123 

MAKES THE DISTANCE FROM THE TARGET TO THE SURFACE A KIMKUM. 

125 

PS=1Q.»^30 126 

XJ = 0.0 127 

YJ = 0.0 128 

ZJ = 0.0 129 

00 16 J=1,NI 12G 

V = ROOT(J) IM 

X2= (XC’A»*2) / ( A»»2-V) 132 

Y2=( YC»A»»2) / (A^»2-V> 13 3 

Z2= (70»B»»2> / (E*^2- V> 13*4 

XA=(-XO»V)/ (A’^2-V) 135 

YA=(-YC»\/) / (A ^^2-V> 136 

ZA= (-70»V)/(R»»2-V> 137 

U=XA»*2+YA»*2+ZA**2 138 

U=U»».5 139 

IF (U.GT.PS) GO TO 16 1 *4 C 

PS=U 1*«1 

XJ = X2 1*42 

YJ=V2 143 

7J = 72 1*4*4 

16 CONTINUF 145 

X2=XJ 146 

Y2 = YJ 1*4? 

7?=ZJ 148 

149 

the altitude equals the MIMMUK distance FFCM THE 150 

TARGET TO THE SURFACE. 151 

CCN=18fl./GK 152 

153 

DETERMINE THE LONGITUOE. 154 

155 

AL0NG=ATAN2( YOfXO) 156 

157 

DETERMINE THE GEODETIC LATITUCE. 158 

159 

D=X2^*2+Y2*»2 160 

D=D'^».5 161 

AA=D»(1.-CE) 162 

ALAT=ATAN2(Z2, AA) 163 

ALAT = SIGN tALAT ,ZO) 164 

Xl=XO*ZP 165 

Yl=YO»ZP 166 

71=Z0*ZP 167 

EX=1./(1.-CE^(SIN(ALAT) »»2)l 168 

EX=A9S(EX) 169 

EX=EX*»,5 170 

EX=EX*A1 171 
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APPENDIX - Continued 



C 



172 


c 


COMPUTE ALTITUDE BY AN EXACT ECUATIOK WHICH FIMMZES PGLND OFF 

173 


c 


ERRORS. 

1 74 

175 

c 



175 




IF(ABS(SIN(ALATJ) .GT..1IGO TO 1<* 

176 




IF(A9S (COS(ALCKG)).GT..l)GO TC 15 

177 




AIT= (Yl/ CCOS (ALAT)*SIN (ALONG) ) )-EX 

176 




GC TO 710 

179 

160 



ALT=(71/SIN( ALAT) )- (EX»(1-CE) ) 

160 




GC TO 710 

181 



15 

ALT = (X1/ (COS( AL AT) *COS( ALONG) ))-EX 

182 



710 

CONTINUE 

133 




ALAT=ALAT»CON 

1 04 

185 

c 



185 


c 


CHANGE SIGN OF LCNGITUCE TO CCNFCFM WITH INPUT WFICH 

106 


c 


WAS IN DEGREES W EST • . . C AtCUL A TI CN S GIVE LONGITLDE IN 

187 


c 


DEGREES EAST, CHANGE SIGN TO GET LONGITUDE IN CEGREES 

188 


c 


WEST, CONVERT FRCH RADIANS TO CEGREES. 

189 

190 

c 



190 




ALONG = -ALONG*CCN 

191 


c 



192 




RETURN 

193 




END 

194 
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APPENDIX “ Continued 


Subroutine QUART 


SUBROUTINE QUART (R, ROO T, NX) 195 

C 196 

C THIS ROUTINE SOLVES FOR THE REAL ZEROS ONLY VIA 197 

C FERRARI'S METHCD, (NEH FIRST COURSE IN THE THEORY OF 198 

5 C EOUATICNSt OICKSCN, PP. 51-52) 199 

C 200 

C ?C1 

C 202 

C R IS THE ARRAY WHICH CONTAINS THE COEFFICIENTS OF THE 203 

10 C QUARTIC arranged IN ASCENOTNG CROER, 20A 

C R(5) * + R(4) * X»*3 ♦ 2C5 

C R(3) * ♦ R(2) » X 4^ F(l) = 0 20E 

C 207 

C THE ARRAY ROOT WILL CONTAIN THE REAL ZEROS CF THE QUARTIC 208 

15 C 209 

C NI IS THE SCALAR VARIABLE THAT STATES THE NUMBER OF REAL ZEROS 210 

C 211 

C 212 

C 217 

20 DIMENSION R(l) ,RCOT(l) 214 

NI=Q 215 

C 216 

C NORMALIZE CCEFFICIENTS SOI 217 

C X»»4 ♦ ♦ C^X»»2 ♦ D»> ♦ E = 0 218 

25 C 219 

OC 19 J=1 f5 220 

19 R(J)=R(J)/R(5) 221 

B=pr4) 222 

C=R(3) 223 

30 0=R(2) 224 

E=R(1) 225 

C 226 

C CALCULATE CCEFFICIENTS fU» Cl* Dl, OF THE RESCLVENT CUBIC 227 

C ♦ B1*Y»»2 4- C1*Y ♦ Cl = 0 228 

35 Bl=-C 229 

C1=B»0-4*E 270 

Q1 = -(B*»2) ♦E4'4^C»E-D»»2 2 31 

C 232 

C IN SOLVING CUBIC EQUATION, CALCULATE COEFFICIENTS 237 

4C C OF THE CORRESPONDING REOUCEC CUBIC WHICH HAS NO TERM 234 

C OF THE SECOND DEGREE BY SETTING Y = Z - 81/3! 235 

C Z»*3 f P*Z 4-0 = 0 236 

C 2^37 

P=Cl-((Rl*»2)/3) 238 

45 Q=01-( (B1»C1) /3)+((2»(01**3)>/27) 239 

C 240 

C THE DISCRIMINANT OF THE GENERAL CUBIC EQUATION IS 2tl 

C EQUAL TO THE DISCRIMINANT CEL CF THE CORRESFONOING 242 

C REDUCED EQUATION. 243 

50 C 244 

DELT = 18»B1»C1*01-4*(B1 »^3) »0l4- (Bl^Cl) ^^2 245 

DELV=-4»(C1*»3) -27*(D1^*2) 246 

0EL=0ELT4-DELV 247 

C 248 

55 C IF OEL IS NEGATIVE* ONE ROOT IS REAL* AND TWO ARE 249 

C CONJUGATE IMAGINARIES, 250 

C 251 


13 



APPENDIX - Continued 





IF(OEL.GE,0, )GC TO 12 

252 




RS=(P/3)^*3+(Q/2)»*2 

253 

60 



RS=ABS(RS) 

254 




A=< (-q/2)*RS»*,5) 

255 




W=((-Q/2)-(RS*^.5)> 

256 




S= ABS(A> 

257 




T=ABS(H) 

258 

65 



W1=W 

259 




A1 = A 

2E0 




VV=1. /3. 

261 




A=S*»VV 

2f2 




W = T*^VV 

2E3 

70 



A=A*SIGN (1, ,A1 ) 

264 




W=W»SIGN(1,,W1) 

2E5 


C 



?£6 


c 


Y1 IS THE SINGLE REAL ROCT CF THE REDUCEC CLRIC» 

267 


c 


Y IS THE SINGLE REAL ROOT CF THE GENERAL CUPIC, 

2E8 

75 

c 



2E9 




Y1=A*W 

270 




Y=Yl-(Bl/3) 

271 




GO TO 10 

27? 



12 

C0N7INLE 

273 

80 

c 



?74 


c 


IF DEL IS POSITIVE, THERE ARE THREE DISTINCT REAL 

275 


c 


ROOTS. IF DEL IS ZERO. THERE ARE AT LEAST TWO EQUAL 

2 76 


c 


REAL ROOTS. 

277 


c 


TRIGONCMETRIC SOLUTION IS USED. 

278 

85 

c 



^79 




TN = { -4^'P/3) 

280 




TN = APS (TN) 

281 




TN=TN^». 5 

282 




C0S3A = -.5*0» ( )»* 1.5> 

283 




IF(AeS(C0S3A UGT.l.) C053A3SIGN( 1..COS3A) 

284 




APC=ACOS (COS3A) /3 

285 




COSA=COS(ARC) 

286 




Yl=TN»COSA 

287 




Y=Yl-(Pl/3» 

288 

95 


10 

CCNTINUE 

2 89 


c 



29 0 


c 


SACK TO SOLVING THE QUARTIC, WHERE Y IS SUCH THAT 

291 


c 


A?»X*»2 f P2»X 4- C2 IS THE SQUARE OF A LINEAR FUNCTION, 

292 


c 


M»X 4 N, AND ALSO EQUAL TO (X^»2 4 B»X/2 f Y/2)*»2 

293 

100 

c 



2 94 




A2=.25»(B»»2)-Ci-Y 

295 




B2=.5*B*Y-D 

296 




C2=.25»(Y^*2» -E 

297 


c 



298 

105 



IF(A2.NE.0.)GC TO 111 

299 


c 



300 


c 


IN CASE A2=0, (X*^2 ♦ B*X/2 ♦ Y/2)**2 = B2»X f C2 

301 


c 


SINCE THE FCLYNCHIAL IS A PERFECT SQUARE, B2 = 0 

302 


c 


AND THE QUADRATICS TO BE SOLVEC ARE 1 

303 

110 

c 


X**2 ♦ B^X/2 4 Y/2 - C2»».5 = 0 ANO 

304 


c 


X»»2 4 B*X/2 4 V/2 ♦ C2*^.5 

305 


c 



306 




DELI = (R/2>»*2 - 4»(Y/2 - C2*^.5) 

307 




DEL2 = (9/2)*»2 - 4*(Y/2 ♦ C2**.5) 

308 
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APPENDIX - Continued 


115 DXl = R/2 ^^09 

0X2 = B/2 
GC TO 223 

111 CONTINUE 312 

C 313 

120 C IN CASE A2 IS NCT 0, THE^ ► = A2»».5, K = B?/(2»H). IIU 

C THEN THE QUADRATICS TO BE SOLVED AREt 315 

C X*»2 ♦ m/?. - M)*X + Y/2 - N = 0 AND 316 

C x»*2 ♦ (e/2 ^ n)*Y ♦ Y/2 ♦ N = 0 317 

C 318 

125 A2=ABS(A2) 319 

AM=A2»»,5 320 

AN=92/(2»AM) 3?1 

DX1=.5»3-AM ’22 

DX2=.5»B+AM 323 

130 DELI = (B/2 - AM>»^2 - 4»(Y/2 - AN) 324 

DEL2 = (B/2 *■ AM)»»2 - 4*(Y/2 ♦ AK) 325 

223 CONTINUE 326 

C 327 

C ROOTS OF TWO QUADRATICS ARE THE FOUR ROOTS CF THE QUARTIC, 328 

135 C CNLY THE REAL ROOTS ARE RETURNED. 329 

C 330 

IF(CEL1.LT.O.)GO TO 22 331 

NI = NH-1 332 

R00T(NI) = J-DXl/2) ♦(( DELl^’.EI / 2) 3 33 

mo NI=NI+1 334 

ROOT (NT) = (-DX1/2 )-( ( DEL1*». 5>/2) 375 

22 IF(CEL2.LT.a.)GO TO 25 336 

NI=NI*1 337 

RCOT(NI)=(-DX2/2J+C(DEL2»*.5)/2) 378 

1C»5 M = KH-1 379 

ROOT(NI) =(-0X2/2) -( (DEL2»'f. 5) / 2) 34 0 

25 CCNTINUF 341 

RETURN 3A2 

END 343 
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APPENDIX - Concluded 


Input and Output Data Samples 


INPUT 

POINT COCPDIh^TES WRT PAOflP SITE 


ELEVATIC^ OEG, 90.0000 
AZI CEG. -0.0000 
RANGE FEET .1000 

RAOAP SITE COORCINATES 

GECCETIC LATITUCE 55.0000 
LONGITUDE ??.0000 
ALTITUDE FEET -0.0000 


CUTPUT 


POINT COORCINATES 

GECCETIC LATITUCE 
LONGITLOE 
ALTITUDE FEET 


55.0000 

22.0000 

.0999 


IhPUT 

POINT COORCINATES HRT RADAR SITE 
ELEVATI0^ DEG. 

AZI DEG. 

RANGE FEET 


90.0000 

- 0.0000 

50000.0000 


RADAR SITE COORDINATES 

GEODETIC LATITUCE 
LONGITUDE 
ALTITUDE FEET 


4^.0000 

121.0000 

500.0000 


OUTPUT- — - 


POINT COORCINATES 

GECCETIC LATITUCE 
LONGITUDE 
altitude FEET 


44.0000 

121.0000 

50499.9900 


INPUT-- — - 

POINT COORCINATES HRT RADAR SITE 
ELEVATION DEG. 

AZI DEG. 

RANGE FEET 


9C.000O 

-O.ODOO 

CCCCCOC99.0000 


RADAR SITE COORDINATES 

GEODETIC LATITUCE 
LONGITUDE 
ALTITUDE FEET 


33.0000 

68.0000 
700.0000 


OUTPUT- — - 


POINT COORCINATES 

GECOETIC LATITUCE 
LONGITUDE 
altitude FEET 


33.0000 

86.0000 

ICG00C0699.O492 


.55C000COn0G6F+02 

.2195«99999fl8E*02 

.998557806015E-DI 


.44COCfl0002€4E+02 
. 120999999999E^03 
. 504999B99952E«'05 


. 329999999979E+02 
. a799999999«8E^02 
. 1COCC0069905E+1D 
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